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ABSTRACT 

We present measurements of carbon, oxygen, silicon, and iron in quasar absorption systems existing 
when the universe was roughly one billion years old. We measure column densities in nine low- 
ionization systems at 4.7 < z < 6.3 using Keck, Magellan, and VLT optical and near-infrared spectra 
with moderate to high resolution. The column density ratios among C II, O I, Si II, and Fe II 
are nearly identical to sub-DLAs and metal-poor ([M/H] < —1) DLAs at lower redshifts, with no 
significant evolution over 2 < z < 6. The estimated intrinsic scatter in the ratio of any two elements 
is also small, with a typical r.m.s. deviation of < 0.1 dex. These facts suggest that dust depletion 
and ionization effects arc minimal in our z > 4.7 systems, as in the lower-redshift DLAs, and that 
the column density ratios are close to the intrinsic relative element abundances. The abundances 
in our z > 4.7 systems are therefore likely to represent the typical integrated yields from stellar 
populations within the first gigayear of cosmic history. Due to the time limit imposed by the age of 
the universe at these redshifts, our measurements thus place direct constraints on the metal production 
of massive stars, including iron yields of prompt supernovae. The lack of redshift evolution further 
suggests that the metal inventories of most metal-poor absorption systems at z > 2 arc also dominated 
by massive stars, with minimal contributions from delayed Type la supernovae or AGB winds. The 
relative abundances in our systems broadly agree with those in very metal-poor, non-carbon-enhanced 
Galactic halo stars. This is consistent with the picture in which present-day metal-poor stars were 
potentially formed as early as one billion years after the Big Bang. 

Subject headings: cosmology: observations — cosmology: early universe — intergalactic medium — 
quasars: absorption lines — stars: abundances 



1. INTRODUCTION 

Metal abundances provide some of our most valuable 
insights into star formation and stellar evolution. The 
elements produced by a single star result from a variety 
of nuclear reactions occurring over its life and death, the 
details of which will depend on the star's mass and ini- 
tial chemical composition, and whether it is a member 
of a binary or higher multiple system. Measurements of 
element abundance ratios can therefore be used to place 
constraints on stellar evolution models. In turn, under- 
standing stellar evolution and nucleosynthetic yields can 
be used to infer the star formation history and mass as- 
sembly of galaxies such as the Milky Way (for reviews, 
see McWilliam 1997; Matteucci 2008). 

Of particular interest are the yields from the first few 
generations of stars. The integrated metal production of 
the so-called "Pop III" stars, which formed out of pristine 
material, will depend strongly on their initial mass func- 
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tion (IMF). In the conventional picture, where Pop III 
stars are predominantly massive due to a lack of metal- 
line cooling in the collapsing gas clouds (e.g., Bromm 
& Larson 2004), their nucelosynthetic signature may be 
highly distinct from later Pop II stars (e.g., Heger & 
Woosley 2002; Nomoto et al. 2006). Recent simulations, 
however, suggest that fragmentation due to turbulence 
may produce a Pop III IMF that is less top heavy than 
previously thought (Clark et al. 2011). In either case, 
the metal abundances in low-mass Pop II stars that sur- 
vive from this period should encode information on the 
nature of the first stars. In addition, early abundance 
patterns should reflect some of the key, and as yet un- 
certain, characteristics governing star formation in the 
early universe, including the IMF. 

A classic approach to studying early star formation is 
to measure metal abundances in very metal-poor stars. 
Significant numbers of objects with mctallicities as low as 

[Fe/H] 5.5 (Christlieb et al. 2002; Frebel et al. 2005) 

have now been identified in the Milky Way and its satel- 
lite galaxies. Indeed, metal-poor stars have produced a 
diverse array of insights, from the nucleosynthetic yields 
of low metallicity stars to the mass assembly of the Milky 
Way (for reviews, sec Beers & Christlieb 2005; Frebel 
2010). 

Using metal-poor stars to determine early abundance 
patterns involves a number of challenges, however. First, 
while highly precise differential stellar abundances can be 
determined under the simplifying assumptions of local 
thermodynamic equilibrium (LTE) and one-dimensional 
model atmospheres, measuring absolute abundances gen- 
erally requires that non-LTE and 3D effects be consid- 
ered. These can introduce large and sometimes uncertain 
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corrections, depending on the element and the choice of 
abundance indicator (for a review, see Asplund 2005). 
Mixing of material within the stellar interior, moreover, 
may contaminate photospheric abundances such that 
they do not reflect the composition of the gas out of 
which the star was formed. Second, although metal-poor 
stars are presumed to be old, the fact that these stars 
tend to be found in the field, rather than in clusters, 
means that age estimates cannot be derived from single 
stellar population modeling. In some cases, radioactive 
dating has confirmed the old ages for metal-poor stars 
with enhanced r-process elements (Hill et al. 2002; Sne- 
den et al. 2003; Frebel et al. 2007a), but this has so far 
been done for only a handful of stars. Uncertainties in 
the ages of most individual metal-poor stars therefore 
remain considerable, although isochrome-bascd methods 
suggest a typical age of ~10-12 Gyr for halo stars with 
-3.0 < [Fc/H] < -1.3 (Jofre & Weiss 2011). Finally, a 
number of metal-poor stars, including a large fraction of 
those with [Fe/H] < —3, display "anomalous" abundance 
features such as strong enhancements in carbon, oxygen, 
and nitrogen (Beers et al. 1992; Beers & Christlieb 2005; 
Cohen et al. 2005; Frebel et al. 2005, 2006; Norris et al. 
2007). The enhanced fraction increases towards lower 
metallicities, with three out of the four known stars with 
[Fe/H] < —4.5 being carbon-enhanced (an exception hav- 
ing only recently been identified by Cafiau et al. 2011). 
The origin of this enhancement is currently debated (e.g., 
Beers & Christlieb 2005; Iwamoto et al. 2005; Meynet 
et al. 2006; Carollo et al. 2011). If it results from mass 
transfer from a binary companion, however (Suda et al. 
2004; Campbell et al. 2010), or if the gas needs to reach 
a threshold level of carbon and/or oxygen before it can 
cool and fragment (Bromm & Loeb 2003; Frebel et al. 
2007b), then the abundances in these stars may not re- 
flect typical ISM abundances at low metallicities. 

Quasar metal absorption lines provide an attractive 
complement to metal-poor stars for measuring primi- 
tive abundances, for several reasons. Quasar absorp- 
tion lines should integrate over larger spatial areas than 
the birthplaces of individual stars, and so may more ac- 
curately reflect mean ISM abundances. A potentially 
more significant advantage, however, is that the conver- 
sion from observed line properties to total element abun- 
dances is often straightforward in quasar absorption sys- 
tems. This is particularly true for the high column den- 
sity absorbers known as damped Lya systems (DLAs, 
logiV(HI) > 20.3), which probe the interstellar media 
of intervening galaxies along quasar lines of sight (e.g., 
Wolfe et al. 2005). These systems are highly optically 
thick to hydrogen-ionizing photons, which ensures that 
nearly all of the elements in the "neutral" -phase gas are 
in the lowest ionization state that has an ionization po- 
tential > 1 Ryd (e.g., carbon, silicon, magnesium, and 
iron will be singly ionized, while oxygen will be fully 
neutral). Furthermore, the densities are low enough that 
nearly all of the ions will be in the ground state. The 
net result is that, in the simplest case, abundances for 
a given element can be computed directly from the col- 
umn density measurements of ions in a single electronic 
state. 7 Complications do arise at lower H I column densi- 

7 For this reason we adopt the common convention among quasar 
absorption line studies of referring to an ion by its absorption spec- 



ties, where ionization effects start to become significant, 
or if there is significant depletion onto dust grains. Lines 
may also become saturated, preventing accurate column 
density measurements. We will return to these points 
later in the paper. 

An equally important advantage of quasar absorption 
lines is the direct access they provide to abundances in 
the early universe. With the discovery of metal-enriched 
absorbers out to redshift six (Becker et al. 2006, 2009, 
2011b; Ryan- Weber et al. 2006, 2009; Simcoc 2006; Sim- 
coe et al. 2011), it is now possible to measure abundances 
in gas within one billion years of the Big Bang. At this 
point, the finite age of the universe requires that the met- 
als we observe must have been produced only by massive 
stars that exploded as core-collapse (Type II and Ib/c) 
and/or prompt Type la supernovae. 8 Indeed, although 
the universe is ~1 Gyr old at z ~ 5-6, current estimates 
suggest that the cosmic stellar mass density at these red- 
shifts doubles on timescales of ^300 Myr (Gonzalez et al. 
2011). It is likely, therefore, that the metals seen in ab- 
sorption at z ~ 5-6 are often less than a few hundred mil- 
lion years old. High-redshift quasar absorption lines thus 
provide unique constraints on nucleosynthesis in massive 
(and presumably metal-poor) stars. At the same time, 
they provide insights into the types of stars that ended 
the cosmic Dark Ages. 

In this paper we present measurements of the relative 
abundances of a-elements (carbon, oxygen, and silicon) 
and iron in a sample of nine quasar absorption line sys- 
tems spanning 4.7 < z < 6.3. These are the highest- 
rcdshift systems for which such measurements have been 
made, and the first within one billion years of the Big 
Bang. They also represent nearly all of the currently 
known low-ionization absorbers at z ~ 5-6, and thus con- 
stitute an unbiased sample with regards to their metal- 
licity. In Section 2 we describe the optical and near- 
infrared spectra used for this work. In Section 3 we 
present the column density measurements for individual 
systems. We then use these to compute relative abun- 
dances in Section 4, and compare the results to lower- 
rcdshift quasar absorption systems. In Section 5 we dis- 
cuss the evidence that the relative abundances in these 
systems represent the typical integrated yields of low 
metallicity, massive stars in the early universe. We also 
compare our measurements to the relative abundances 
in metal-poor Galactic halo stars. We summarize our 
conclusions in Section 6. 

Throughout the paper we express logarithmic relative 
abundances with respect to the solar values as [X/Y] = 
log (Nx/Ny) — log(nx/^Y)o, where Nx is the column 
density of element X. All quantities are computed with 
respect to the solar photospheric values from Asplund 
et al. (2009). 



2. THE DATA 

trum. Hence, N(C II) is the column density of carbon ions pro- 
ducing C II absorption, or cquivalcntly, the column density of C+ 
ions. 

8 Here we distinguish between a potential "prompt" contribution 
to the Type la rate from high-mass stars, and a "delayed" contri- 
bution from evolved low-mass stars (Dallaporta 1973; Oemler & 
Tinsley 1979; Mannucci et al. 2005; Scannapieco & Bildsten 2005; 
Matteucci et al. 2006). 
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TABLE 1 
Data 



QSO 




-cm 


Instrument 


^cxp 










(hrs) 


SDSS J0040- 


-0915 


4.98 


MIKE 


8.3 








NIRSPEC 


11.0 


SDSS J0231- 


-0728 


5.42 


X-Shooter 


6.0 


SDSS J1208- 


-0010 


5.27 


X-Shooter 


9.0 


SDSS J0818- 


-1722 


6.00 


NIRSPEC 


25.5 


SDSS J1148- 


-5251 


6.42 


NIRSPEC 


7.5 



Note. Keck/HIRES optical spectra for 

SDSS J0231-0728, SDSS J0818+1722, and 
SDSS J1148+5251 were presented in Becker et al. 
(2006, 2011b). 



The data used for this work include optical and 
near-infrared spectra of five QSOs at 4.98 < z cm < 
6.42. Metal absorption line measurements from 
Keck/HIRES (FWHM = 6.7 km s" 1 ) optical spec- 
tra of SDSS J0231-0728, SDSS J0818+1722, and 
SDSS J1148+5251 were presented in Becker et al. (2006, 
2011b). The remaining observations are summarized in 
Table 1. Optical data from Magellan/MIKE (FWHM 
= 13.6 km s- 1 ) for SDSS J0040-0915 were also used in 
Becker et al. (2011a); however, the metal absorption lines 
are presented for the first time here. 

We obtained near-infrared spectra of 
SDSS J0040-0915, SDSS J0818+1722, and 
SDSS J1148+5251 with Keck/NIRSPEC (McLean 
et al. 1998) between December 2009 and March 2010. 
Each object was observed using a single setting in 
either the Nirspec-5 (1.41-1.81 /im; SDSS J0040-0915 
and SDSS J0818+1722) or Nirspec-6 (1.56-2.32 fj,m; 
SDSS J1148+5251) band. The settings were chosen 
to cover at least one strong Fe II line. The data were 
reduced using a suite of custom routines, described 
in Becker et al. (2009). These were designed to op- 
timally subtract the sky background and extract the 
one-dimensional spectra, while minimizing the noise 
due to a variety of detector issues. Corrections for 
telluric absorption were made using the atmospheric 
transmission spectrum of Hinkle et al. (2003). The 
spectral resolution was measured from telluric absorp- 
tion lines in standard star spectra, and was found to be 
-15-16 km s" 1 (FWHM). 

Optical and near- infrared spectra of SDSS J023 1—0728 
and SDSS J1208+0010 were obtained with VLT/X- 
Shooter (D'Odorico et al. 2006) between October 2009 
and March 2010. Slit widths of 0'.'9 were used in both 
the VIS and NIR arms. The data were reduced using a 
custom pipeline similar to that used for the NIRSPEC 
data, and atmospheric absorption corrections were again 
made using a model. The resolution, measured from tel- 
luric lines in the quasar spectra, was found to be ~23 
(~32) km s- 1 in the VIS arm, and -31 (-34) km s" 1 in 
the NIR arm for SDSS J0231-0728 (SDSS J1208+0010). 
Resolutions higher than the nominal values were achieved 
due to good seeing. 

3. COLUMN DENSITY MEASUREMENTS 

Column densities for C II, O I, and Si II were mea- 
sured from optical spectra, while Fe II measurements 
were mostly made from strong lines in the infrared. The 



measurements were done using two techniques. In most 
cases where the system appeared to have a single, narrow 
component, Voigt profiles were fit using vpfit. When 
fitting, we explicitly tied together the Dopplcr parame- 
ters for all ions, which assumes the lines are turbulently 
broadened. For the remaining systems we used the ap- 
parent optical (AOD) method (Savage & Sembach 1991), 
which calculates a total column density by summing the 
optical depth over all pixels within a specified veloc- 
ity interval. For noisy data (e.g., the HIRES data for 
SDSS J0231-0728), we first smoothed the spectrum by 
convolving the flux with a Gaussian kernel with FWHM 
equal to one half the instrumental resolution before mea- 
suring optical depths. This has only a minor effect on the 
spectral resolution, but prevents individual pixels from 
dominating the optical depth by taking advantage of the 
fact that the absorption lines are expected to have a finite 
width. For further discussion, see Becker et al. (2011b). 

Notes on individual systems are given below. A sum- 
mary of the measured column densities is given in Ta- 
ble 4. 

3.1. SDSS J0040-0915, z = 4.7393 

This system shows mildly saturated C II A1334 and 
O IA1302, from which we can only derive lower limits on 
the column densities. Multiple weak or moderate lines 
of Si II and Fe II are present, however. We also cover 
Mg II in the NIRSPEC data, although it is saturated. 
Moderate C IV and Si IV absorption appears, but these 
lines are likely to arise from a separate, high-ionization 
phase (e.g., Fox et al. 2007). The metal lines for this 
system are shown in Figure 1. At this redshift, the Lya 
forest is highly absorbed. We can, however, set an upper 
limit on the H I column density for this system from Lya 
and Ly/3 (log V H i < 20.4 (cm" 2 ); Figure 2). 

3.2. SDSS J1208+0010, z = 5.0817 

This system appears to have a single component. It 
is unresolved in the X-Shooter data, however, and we 
lack supplementary high-resolution spectra. We have 
therefore attempted to obtain column densities by fit- 
ting Voigt profiles, with the 6-parameter constrained by 
simultaneously fitting multiple Si II lines with differ- 
ent oscillator strengths. A single-component fit gives 
b = 9.5 ± 0.2 km s _1 , with which column densities for 
all ions can then be derived. The column densities for 
Si II and Fe II are likely to be robust, as we have at least 
one line in each case that are on the linear part of the 
curve of growth. C II and O I are more problematic, 
however, as there could be hidden saturation issues. We 
found that we could match the observed profiles of these 
ions with multiple, blended components that had a to- 
tal column density significantly greater than that of the 
best-fit single-component fit. Lower limits are therefore 
given for C II and O I, as well as for Mg II, for which 
the A2796 line is strongly affected by sky-line residuals. 
Weak high-ionization lines (C IV and Si IV) are detected 
for this system. All of the metal lines are plotted in Fig- 
ure 3. The regions of the forest covering Lya and Ly/3 are 
shown in Figure 3. The upper limit on the H I column 
density is log N H , < 20.3 (cm" 2 ). 

3.3. SDSS J02S1-0728, z = 5.3380 
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Fig. 1. — Metal lines for the z = 4.7393 system towards SDSS J0040— 0915. Data covering transitions with rest wavelength up to 1608 A 
are from MIKE, while the remainder are from NIRSPEC. The normalized flux for each ion is plotted at the nominal rcdshift of the system. 
Shaded regions show detected transitions. Peaks in the error array, which is plotted along the bottom of each panel, indicate regions 
affected by skylines, for example in the red wings of Si II 1526 and C IV 1548. 
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Fig. 2. — Lya and Ly/3 transitions for the z = 4.7393 system 
towards SDSS J0040-0915. The data are from MIKE. Absorption 
profiles are not fits, but show the maximum allowed H I column 
density, log JV H 1 = 20 .4 (cm -2 ). 



The HIRES data for this system was originally pre- 
sented in Becker et al. (2006). We have re- reduced the 
data to take advantage of a number of improvements in 
the pipeline, and have re-measured column densities for 
C II, O I, and Si II. Lines of Fe II and Mg II, as well 
Si II A 1526, are covered in the X-Shooter data. Fe II is 
detected in two lines with high oscillator strengths, A2344 
and A2382. We note that these lines are often badly sat- 



urated in lowcr-redshift DLAs and sub-DLAs, but are 
quite weak here. Extended C IV and Si IV are present, 
but these are weak at the velocity of the strongest low- 
ionization component, where O I is detected. The metal 
lines are plotted in Figure 5. The upper limit on the H I 



column density is logiVni < 20.3 
2006). 



(cm 2 ) (Becker et al. 



3.4. SDSS J0818+1722, z = 5.7911 

The HIRES data for this system were presented in 
Becker et al. (2011b). Fe II is detected in the NIR- 
SPEC data for two lines, A2344 and A2382. Although 
the absorption is relatively weak, these show the same 
extended velocity structure present in the other low- 
ionization lines (Becker et al. 2011b). The Fe II lines 
are shown in Figure 6. 

3.5. SDSS J 08 18+ 1722, z = 5.8765 

The HIRES data for this system were presented in 
Becker et al. (2011b). Although we cover Fe II A2344 
and A2382 in the NIRSPEC data, no significant absorp- 
tion is detected. The expected strongest line, A2382, un- 
fortunately falls in a region of the spectrum with a low 
signal-to-noise ratio. An upper limit is therefore given 
for logTVpeo in Table 4. The regions of the spectrum 
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Fig. 3. — Metal lines for the z = 5.0817 system towards SDSS J1208+0010. All data are from X-Shooter. Lines are as in Figure 1. 
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Fig. 4. — Lyct and Ly/3 transitions for the z = 5.0817 system 
towards SDSS J1208+0010. The data arc fromX-Shooter. Ab- 
sorption profiles are not fits, but show the maximum allowed H I 
column density, logTVni = 20.3 (cm~ 2 ). Although the flux ex- 
ceeds the absorption profile far out in the red wing of Lya, this 
discrepancy strongly depends on the uncertain continuum level. 

covering these lines are plotted in Figure 7. 

3.6. SDSS Jl 148+ 5251, z = 6.0115 

The HIRES data for this system were presented in 
Becker et al. (2006) and re-analyzed in Becker et al. 
(2011b). Fe II A2600 is detected in the NIRSPEC data. 
Although it falls adjacent to a region strongly contami- 
nated by skylines, the Fe II line is not affected, as shown 
in Figure 8. 



3.7. SDSS J1148+5251, z = 6.1312 

The HIRES data for this system were presented in 
Becker et al. (2006) and re-analyzed in Becker et al. 
(2011b). We clearly detect Fe II A2382 in the NIRSPEC 
data. This line is plotted in Figure 9. 

3.8. SDSS J1148+525L z = 6.1988 

The HIRES data for this system were presented in 
Becker et al. (2006) and re-analyzed in Becker et al. 
(2011b). This is the weakest of our z > 4.7 systems, and 
although we cover Fe II A2344 in the NIRSPEC data, the 
spectrum is too noisy at the corresponding wavelength to 
allow any useful constraint on the column density. 

3.9. SDSS J1148+ 5251, z = 6.2575 

The HIRES data for this system were presented in 
Becker et al. (2006) and re-analyzed in Becker et al. 
(2011b). We cover Fe II A2382 in the NIRSPEC data, but 
no significant absorption is detected. The corresponding 
region of the spectrum is shown in Figure 10. 

4. RELATIVE ABUNDANCES 

The Lya forest at z > 5 is too highly absorbed to allow 
accurate H I column density measurements. At z > 6, 
the forest becomes almost completely opaque, and H I 
cannot be measured at all. We are therefore unable to 
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Fig. 5. — Metal lines for the z = 5.3380 system towards SDSS J0231— 0728. Data covering transitions with rest wavelength up to 
1334 A are from HIRES, and were re-reduced from Becker et al. (2006). The remaining data are from X-Shooter. Lines are as in Figure 1. 
C IV 1551 is blended with unrelated absorption lines. 



TABLE 2 

Relative Abundances of Individual z > 4.7 Systems 



QSO 


^abs 


[Si/O] 


[C/O] 


[C/Si] 


[C/Fc] 


[O/Fe] 


[Si/Fe] 


[Mg/Fe] 


SDSS J0040-0915 


4.7393 


< 0.33 




> -0.46 


> -0.12 


> 0.02 


0.35 ± 0.06 


> -0.45 


SDSS J1208+0010 


5.0817 


< 0.27 




> -0.33 


> 0.14 


> 0.20 


0.47 ± 0.07 


0.28 ±0.13 


SDSS J0231-0728 


5.3380 


-0.14 ±0.06 


-0.42 ± 0.07 


-0.29 ±0.06 


0.13 ±0.06 


0.55 ±0.06 


0.42 ± 0.05 


0.20 ± 0.04 


SDSS J0818+1722 


5.7911 


0.00 ± 0.05 


-0.15 ±0.04 


-0.15 ±0.05 


0.30 ±0.07 


0.45 ±0.07 


0.45 ± 0.08 




SDSS J0818+1722 


5.8765 


-0.08 ±0.07 


-0.17 ±0.09 


-0.09 ±0.09 


> -0.02 


> 0.15 


> 0.07 




SDSS J1148+5251 


6.0115 


0.04 ± 0.04 


-0.25 ±0.06 


-0.29 ±0.06 


-0.00 ±0.25 


0.25 ±0.25 


0.29 ± 0.25 




SDSS J1148+5251 


6.1312 


-0.32 ±0.24 


-0.65 ±0.27 


-0.34 ±0.16 


-0.09 ±0.23 


0.56 ±0.29 


0.25 ± 0.19 




SDSS J1148+5251 


6.1988 


-0.20 ±0.14 


-0.34 ±0.17 


-0.14 ±0.12 










SDSS J1148+5251 


6.2575 


-0.02 ± 0.17 


-0.08 ± 0.26 


-0.06 ±0.23 


> 0.13 


> 0.21 


> 0.18 





Note. — Relative abundances are calculated as [X/Y] = log (N-^/Ny) — log (nx/riY)©> without any corrections for dust or ionization 
effects. The column densities used are those integrated over the velocity interval spanned by O I (see Table 4). Solar values are from 
Asplund ct al. (2009). 



directly estimate metallicities for our systems. We can, 
however, estimate element ratios based solely on the ionic 
column densities. 

Directly converting from an ionic column density to a 
total element column density (log iVsi n to log Nsi, for ex- 
ample) requires first that there is minimal depletion onto 
dust grains, and second, that the ionization corrections 
are known. Dust-depletion in DLAs tends to increase 
with metallicity, but is negligible at [Fe/H] < — 2 (Vladilo 
2004). The trend of declining DLA metallicity with red- 
shift (e.g., Wolfe et al. 2005) suggests that our z > 4.7 
systems may be in this very metal poor regime. Never- 
theless, differential dust depletion is a potential source 
of error in determining relative abundances, an issue to 



which we will return below. 

Ionization corrections are a potentially greater cause 
for concern. In the simplest case, the corrections will be 
small and can be safely neglected. This is true for heavily 
self-shielded systems, such as DLAs, in which nearly all 
of the hydrogen is in H I, and nearly all of the metals are 
in their neutral or first excited states. Oxygen has a first 
ionization potential very similar to that of hydrogen, and 
would therefore expected to be present as O I. Silicon, 
carbon, iron, and magnesium, meanwhile, have first ion- 
ization potentials below that of hydrogen, and will not 
be shielded by neutral gas. These should therefore be 
singly ionized (for discussions of ionization corrections 
in DLAs see, e.g., Prochaska et al. 2002; Pettini et al. 
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Fig. 6.— NIRSPEC data covering Fc II for the z = 5.7911 
system towards SDSS J0818+1722. Lines are as in Figure 1. Data 
covering other ions for this system were presented in Becker et al. 
(2011b). 
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Fig. 7.— NIRSPEC data covering Fe II for the z = 5.8765 system 
towards SDSS J0818+1722. The absorption for Fe II A2382 is not 
a significant detection. Lines are as in Figure 1. Data covering 
other ions for this system were presented in Becker et al. (2011b). 
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Fig. 8.— NIRSPEC data covering Fe II for the z = 6.0115 system 
towards SDSS J1148+5251. The Fe II 2600 line falls redward of a 
region strongly affected by skylines, and is clearly detected. Lines 
are as in Figure 1. Data covering other ions for this system were 
presented in Becker et al. (2006, 2011b). 
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Fig. 9.— NIRSPEC data covering Fe II for the z = 6.1312 
system towards SDSS J1148+5251. Lines are as in Figure 1. Data 
covering other ions for this system were presented in Becker et al. 
(2006, 2011b). 

2008). For sub-DLAs (19.0 < logA Hl < 20.3), the H I 
self-shielding is significantly reduced, and ionization cor- 
rections up to ^0.3 dex for an individual ion have been 
estimated (e.g., Peroux et al. 2007). An important ques- 
tion, therefore, is whether significant ionization effects 
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Fig. 10.— NIRSPEC data covering Fe II for the z = 6.2575 
system towards SDSS J1148+5251. The absorption for Fe II 2383 is 
not a significant detection. Lines are as in Figure 1. Data covering 
other ions for this system were presented in Becker et al. (2011b). 

may be present in our z > 4.7 systems. We unfortunately 
do not have useful constraints on any intermediate ion- 
ization states (e.g., Si III, Fe III, or Al III) because these 
transitions are either lost in the forest, are not covered 
by our data, or are too weak to be detected. We note, 
however, that when inferring relative abundances, it is 
the difference in the ionization correction between two 
elements that matters. Although different elements will 
be ionized at different rates, the effects will almost always 
act in the same sense, which is to reduce the population 
of the "low" (O I, Si II, C II, etc.) ionization state. Even 
for a partially ionized system, therefore, the ratio of two 
ions (Si II and Fe II, for example) may be expected to be 
reasonably close to the total element ratio. 

For this work, we will use the fact that all means of 
affecting ionic column density ratios-dust depletion, ion- 
ization effects, as well as intrinsic variations in the rel- 
ative abundances-will then tend to increase the scatter 
between systems. The presence of O I, which has an ion- 
ization potential close to 1 Ryd, already signals that ion- 
ization corrections for our z > 4.7 systems may be small. 
We will argue, moreover, that the low level of scatter 
in the column density ratios, and their agreements with 
the ratios in lower-rcdshift, metal-poor DLAs, strongly 
suggests that they reflect the underlying relative abun- 
dances. 

The uncorrected relative abundances for individual 
z > 4.7 systems are given in Table 2. These do not 
include any corrections for dust depletion or ionization 
effects. We compare the z > 4.7 results to values for 
lower-redshift systems in Figure 11. The lower-redshift 
data include three groups of low-ionization absorption 
systems measured from high-resolution data, for which 
accurate column density measurements for some of the 
most abundant ions (O I, Si II, C II, and Fe II) are avail- 
able. The first group is a compilation of very metal- 
poor DLAs from Cooke et al. (2011b), which includes 
their own measurements and other results from the lit- 
erature (see references therein). These DLAs were se- 
lected to have very low mctallicities ([Fe/H] < —2.0), 
and hence weak absorption lines. In particular, column 
densities have been measured for O I and C II, whose 
transitions arc normally saturated for higher mctallicity 
DLAs (e.g., Prochaska k Wolfe 2002). We omit the po- 
tentially carbon-enhanced metal-poor DLA from Cooke 
et al. (2011a) due to possible large uncertainties in the 
column densities. The velocity width of this system is 
sufficiently narrow (b = 2.36 km s _1 ) that thermal broad- 
ening may be significant. The fit is largely driven by the 
multiple lines of Si II, for which the corresponding tem- 
perature would be ~9500 K. In a thermally broadened 
model, the C II column density would be substantially 
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Fig. 11. — Uncorrected relative abundances in low-ionization systems as a function of redshift. In all cases, the relative abundances are 
determined directly from the column density measurements, and are not corrected for dust depletion or ionization effects. Filled circles 
are from this work. Open circles are sub-DLAs from Dessaugcs-Zavadsky et al. (2003) and Pcroux ct al. (2007). Triangles are metal-poor 
DLAs from Cooke et al. (2011b), and references therein. Squares are DLAs with [M/H] < —1.0 from Wolfe et al. (2005). Arrows are 1-a 
upper or lower limits, which at z > 4.7 are for this work and otherwise are for literature values. The mean value with a 1-tr uncertainty is 
given in each panel, where the values have been calculated from all measurements (see Table 3). Estimates of the nominal intrinsic r.m.s. 
scatter, <rn> an d the 95% upper limit on <tq (in parentheses) arc shown as dark and light shaded bands, respectively. 



(> 1 dex) lower than in the turbulent model they as- 
sume (confirmed by R. Cooke, private communication). 
We note that thermal broadening is much less likely to 
be an issue for our systems, as our narrowest absorber 
has b = 3.8 km s™ 1 , for which the corresponding temper- 
ature for Si II would be ^25000 K. We further compare 
our z > 4.7 data to a sample of DLAs from Wolfe et al. 
(2005) with metallicities -2.7 < [Si/H] < -1.0. These 
have column densities measured for Si II and Fe II, which 
have weak transitions that remain optically thin. The 
upper cutoff in metallicity is chosen to avoid systems 
with an enhancement in [Si/Fe] that is most plausibly 
due to dust depletion (Wolfe et al. 2005). Finally, we in- 



clude a selection of sub-DLAs from Dessauges-Zavadsky 
et al. (2003) and Peroux et al. (2007). For these systems 
we restrict ourselves to measurements where the metal 
lines fall redward of the Lya forest, are not strongly sat- 
urated, and are free from any apparent contamination 
from other lines. We have not included the metal-poor 
DLA measurements from Penprase et al. (2010), since 
these were made from medium-resolution spectra. Al- 
though they are broadly consistent with the other lit- 
erature values, the trends in relative abundances with 
metallicity identified by Penprase et al. could constitute 
an intriguing counterexample to the results described be- 
low if confirmed with high-resolution data. 
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TABLE 3 
Mean Relative Abundances 



Group 


# Systems Mean 






[Si/O] 


sub-DLAs 


1 






VMP-DLAs 


19 


-0.08 ±0.02 


0.09 (< 0.13) 


z > 4.7 Systems 


7 


-0.04 ±0.03 


0.05 (< 0.18) 


All 


27 


-0.08 ±0.02 


0.08 (< 0.11) 


[C/O] 


sub-DLAs 


2 






VMP-DLAs 


8 


-0.30 ±0.05 


0.11 (< 0.22) 


2 > 4.7 Systems 


7 


-0.26 ±0.05 


0.09 (< 0.30) 


All 


17 


-0.28 ±0.03 


0.10 (< 0.18) 


[C/Si] 


sub-DLAs 


2 






VMP-DLAs 


8 


-0.18 ±0.02 


0.04 (< 0.12) 


2 > 4.7 Systems 


7 


-0.21 ±0.04 


0.04 (< 0.17) 


All 


17 


-0.19 ±0.02 


0.03 (< 0.08) 


[C/Fe] 


sub-DLAs 


1 






VMP-DLAs 


6 


0.14 ±0.03 


0.04 (< 0.29) 


z > 4.7 Systems 


1 


0.17 ±0.07 


0.09 (< 0.47) 


All 


11 


0.15 ±0.03 


0.04 (< 0.16) 


[O/Fe] 


sub-DLAs 


1 






VMP-DLAs 


12 


0.39 ±0.02 


0.04 (< 0.08) 


z > 4.7 Systems 


1 


0.50 ± 0.05 


0.00 (< 0.36) 


All 


17 


0.40 ± 0.02 


0.04 (< 0.08) 


[Si/Fe] 


sub-DLAs 


9 


0.40 ± 0.04 


0.09 (< 0.24) 


Metal-poor DLAs 


31 


0.38 ± 0.02 


0.10 (< 0.13) 


VMP-DLAs 


17 


0.33 ±0.02 


0.06 (< 0.10) 


z > 4.7 Systems 


6 


0.41 ±0.03 


0.00 (< 0.12) 


All 


63 


0.37 ±0.01 


0.09 (< 0.11) 



Note. — Relative abundances are computed assum- 
ing no dust depletion or ionization corrections. Mean 
values are given only for samples with > 4 systems. 

a Estimated intrinsic scatter in the relative abundances. 
The 95% confidence upper limit on ctq is given in paren- 
theses. 



The mean uncorrected relative abundances for the 
z > 4.7 and lowcr-rcdshift systems are given in Table 3. 
We give the relative values for all six combinations of our 
four measured elements, though these are obviously not 
all independent. Limits are not included, though these 
are generally consistent with the mean values. We also 
estimate the intrinsic scatter in the relative abundances 
by assuming a log-normal distribution and computing 
the r.m.s. logarithmic scatter, a$, which, when added 
in quadrature to the measurement errors, produces a re- 
duced x 2 equal to 1.0 for a fit to a constant value. This 
nominal value of oo is added to the individual errors when 
computing the mean values. A 95% upper limit on <7o 
is also estimated as that which produces the 95% lower 
bound on x 2 for the available number of measurements. 

Two points stand out from Figure 11 and Table 3. 
First, the relative abundances are very similar among 
different groups of absorbers. For a given element ratio, 
all of the mean values agree to within < 2a, and most 
agree to within < la. Second, the intrinsic scatter for all 
element ratios appears to be small. The nominal intrinsic 



scatter within the complete sample is always < 0.1 dex 
(r.m.s.). Even the 95% upper limits on ao are gener- 
ally < 0.1. The exceptions are [C/O], for which both 
C II and O I are measured from only one or two strong 
lines, and [C/Fe], where the number of data points with 
small errors is low. Values of <To larger than 0.1 are also 
generally allowed for individual groups of systems only 
when there are relatively few measurements and/or the 
measurement error bars are large, and so adding intrinsic 
scatter has less impact on the total scatter. 

We emphasize that the intrinsic scatter will include 
any variations in dust depletion and ionization correc- 
tions, and should therefore give an upper limit on the 
true variation in the element ratios. The intrinsic scat- 
ter may also be overestimated if the measurement un- 
certainties are too small. Absorption line measurements 
are susceptible to errors caused by contamination from 
unrelated lines, which may be difficult to detect. This 
is a particular problem when the measurements rely on 
only one or two lines, as is always the case for O I and 
C II, and may also be true for Si II and Fe II, depend- 
ing on the strength of the system and the wavelength 
coverage of the data. Hidden saturation effects, even in 
high-resolution data, arc another potential source of un- 
certainty that may not be properly included in the formal 
measurement errors. These factors all suggest that the 
minimal scatter we measure in the column density ratios 
strongly indicates that the variation in the underlying 
element ratios is genuinely small. The lack of strong 
variations, in turn, suggests that the column density ra- 
tios are a good indication of the true relative element 
abundances. 

5. DISCUSSION 

We have shown that the ionic column density ratios in 
our z > 4.7 systems are remarkably similar to those in 
low-ionization systems at z ~ 2-4. By necessity, most of 
the lower-rcdshift data are drawn from a sample of low- 
metallicity DLAs; these have relatively weak metal lines, 
making it possible to measure column densities for O I 
and C II. While the metal-poor systems are a convenient 
reference sample, it should be remembered that these 
were deliberately drawn from the tail of the metallic- 
ity distribution, and may not represent typical DLAs at 
these redshifts (although typical DLAs with metallicitics 
up to [Si/H] = —1.0 are included in the [Si/Fc] measure- 
ments). In contrast, our z > 4.7 sample includes nearly 
all of the low-ionization systems detected in an unbiased 
survey (Becker et al. 2006, 2011b, and additional obser- 
vations of z > 5 quasars by our group). We turn now to 
consider the possible origins of these high-redshift metals, 
and to compare our results with the relative abundances 
in metal-poor stars. 

5.1. Pop II or Pop III Enrichment? 

It is worth considering whether our high-redshift sys- 
tems reflect the yields of conventional Pop II stars, or 
perhaps an early generation of metal-free stars. Cooke 
et al. (2011b) noted that the relative abundances in their 
metal-poor DLAs were consistent with enrichment from 
either low-metallicity {Z < 1/3 Zq) Pop II stars or zero- 
metallicity "Pop III" stars (with masses of 10 to 100 M Q ), 
based on comparisons with theoretical yields. This sug- 
gests that our systems could also be dominated by met- 
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als from either population. There are several issues to 
note, however. First, much of the discriminating power 
in the lower-redshift systems came from the abundances 
of nitrogen and aluminum, which we are unfortunately 
unable to measure. 9 Equally important, however, are the 
uncertainties in the theoretical yields. These depend on 
a broad array of minimally constrained factors, includ- 
ing mixing, mass loss, and the mechanics of supernovae. 
The yields are also highly sensitive to stellar mass (e.g., 
Woosley & Weaver 1995; Heger & Woosley 2010), which 
means that interpolating over a sparsely sampled mass 
grid (e.g., Chieffi & Limongi 2004; Nomoto et al. 2006) 
may not give accurate mean abundances. In addition, 
the results will depend on the lower and upper mass lim- 
its over which the yields are integrated. 

We suggest that the abundance ratios in our high- 
redshift systems may be most useful as empirical con- 
straints on the integrated yields of (presumably low 
metallicity) massive Pop II stars. Although the first gen- 
eration of stars must necessarily be metal-free, the phase 
of cosmic star formation dominated by Pop III stars is 
expected to be brief, as high-density regions will quickly 
become enriched above the critical metallicity threshold. 
The transition to Pop II stars should occur well before 
z ~ 10 (e.g., Maio et al. 2010). Current estimates of the 
buildup of stellar mass at z > 4 suggest that only ~l/3 
of the stellar mass at z ~ 6 would have been in place by 
z ~ 10 (Gonzalez et al. 2011). It is therefore reasonable 
to assume that Pop II stars have dominated the metal 
enrichment by z ~ 6. As noted above, our high-redshift 
systems represent typical metal-enriched regions of the 
universe at z ~ 5-6, at least when weighting by cross 
section. Their relative abundances are therefore likely to 
reflect the integrated yields of Pop II stars, and massive 
stars in particular, since there has not been enough time 
for low-mass stars to evolve. Over longer timescales, con- 
tributions from lower-mass stars through delayed Type la 
supernova (which will add iron) and AGB winds (which 
will add carbon) will start to become important. The 
similarity in the abundances between our z > 4.7 sys- 
tems and those at 2 < z < 4.5, however, suggests that 
supernovae of massive stars remain the dominant source 
of enrichment even at lower redshifts. This may either be 
because the moderate redshift, low- metallicity systems 
were enriched at earlier times, or because they were en- 
riched relatively recently. 

5.2. Comparison with Metal-Poor Stars 

Finally, we can compare the abundances in the high- 
redshift absorbers to those in metal-poor stars. Figure 12 
shows the ratios of C, O, and Si to Fe as a function of 
metallicity for a sample of very metal-poor halo giant 
stars from Cayrel et al. (2004). For each element, we also 
plot the estimated intrinsic Ter range in [X/Fe] spanned 
by the quasar absorption line systems. The agreement 
is generally good, which is consistent with the picture 
in which present-day metal-poor stars were potentially 
formed at very high redshifts. It should be noted, how- 
ever, that the stellar values depend on the choice of model 

9 The transitions for N I and N II fall in the Lyo forest, which is 
highly absorbed at z > 4.7. While Al II A1670 can in principle be 
measured for our systems, in nearly all cases where we had spectral 
coverage the data were of too low quality to place meaningful limits. 
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Fig. 12. — Ratios of carbon, oxygen, and silicon to iron for the 
sample of metal-poor stars from Cayrel et al. (2004). All stellar 
values have been adjusted to Asplund et al. (2009) solar values. 
For carbon, only "unmixed" stars are shown (Spite et al. 2005). 
Values for oxygen include corrections for three-dimensional effects 
(see notes in Cayrel et al 2004). The carbon-enhanced stars CS 
22949-037 and CS 22892-052 are plotted with open squares. In 
each panel, the values for low-ionization quasar absorption line 
systems are shown as shaded bands. Dark grey bands span the 
nominal l-cr r.m.s. intrinsic scatter (cro) around the mean value 
in the absorption-line data, while light grey bands reflect the 95% 
upper limit on <ro (see Figure 11 and text for details). 

parameters. The stellar carbon abundances were deter- 
mined from the G-band of CH at ^4300 A. These are 
1D/LTE values, although the (negative) corrections due 
to 3D and non-LTE effects are potentially large (Asplund 
2005; Collet et al. 2006; Frebel et al. 2008). Atmospheric 
carbon abundances in giant stars can be altered by con- 
vective mixing. We therefore plot [C/Fe] values only for 
those stars that appear to be unmixed based on their rel- 
ative abundances of C, N, and Li (Spite et al. 2005). Oxy- 
gen was measured from the forbidden [O I] 6300 A line, 
which is not expected to suffer non-LTE effects. Three- 
dimensional effects can be significant at low mctallicitics, 
however, and we plot the Cayrel et al. (2004) values for 
which a ID to 3D correction of ~ —0.3 has been applied, 
following the results of Nisscn et al. (2002). Silicon abun- 
dances were computed from a single line at 4103 A, and 
arc 1D/LTE. 



Fe and a-elements at z ~ 5-6 



11 



The C/0 ratio is of particular interest as a diagnos- 
tic of the IMF, as more massive stars tend to produce 
more oxygen relative to carbon. An added advantage 
is that both carbon and oxygen are formed during the 
hydrostatic burning phase, which helps to simplify the 
comparison with theoretical yields. The most direct con- 
straints on the IMF can potentially be placed at low 
metallicities, where mass loss from stellar winds is mini- 
mal. The results from stellar abundances, however, again 
depend somewhat on the choice of model parameters. 
Tomkin et al. (1992) and Akerman et al. (2004) mea- 
sured [C/O] in metal-poor dwarf and giant stars, respec- 
tively, using high-excitation C I and O I lines. These 
lines have the advantage that, while the separate C and 
O measurements both depend strongly on the effective 
temperature, their ratio generally does not. Non-LTE 
effects were included for the dwarf stars by Tomkin et al. 
(1992) and for the giant stars by Fabbian et al. (2009). 
Three-dimensional modeling was not done in either case, 
but the effects are expected to be small (Fabbian et al. 

2009). Both groups found that [C/O] 0.6 to -0.5 

over -2.0 < [O/H] < -0.5 (-2.5 < [Fc/H] < -1). 
The Akerman et al. (2004) sample extends to some- 
what lower metallicities, and exhibits an upturn at the 
low-metallicity end, reaching [C/O] ~ —0.3 to 0.0 at 
[O/H] ~ —3, depending on how the non-LTE corrections 
are calculated (Fabbian et al. 2009). This trend is also 
seen in the metal poor giants studied by Cayrel et al. 
(2004) and Spite et al. (2005), but only when the nega- 
tive 3D corrections are applied to the oxygen abundances 
(and not to the carbon). 

The apparent upturn in stellar [C/O] values at the 
lowest metallicities has been cited as evidence of early 
enrichment from Pop III stars and/or a top-heavy IMF 
(Akerman et al. 2004). Very metal-poor (Cooke et al. 
2011b) and very high-redshift (this work) quasar absorp- 
tion line systems have [C/O] ~ —0.3, similar to the most 
metal-poor stars. The implications for the IMF, how- 
ever, will depend on the adopted theoretical yields (see 
also the discussion in Cooke et al. 2011b). For example, 
[C/O] ~ —0.3 is close to the value calculated by Cescutti 
et al. (2009) for a standard Kroupa (2002) IMF corrected 
for binaries, using the yields of Meynet & Maeder (2002), 
which include the effects of rotation and mass loss (these 
yields were also used by Akerman et al. 2004). Based on 
this calculation, the [C/O] ratios seen in quasar absorp- 
tion line systems and the most metal-poor stars show 
no evidence for a top-heavy IMF. In that case, how- 
ever, it becomes challenging to explain why even lower 
[C/O] ratios are seen in galactic halo and bulge stars 
with —2.0 < [O/H] < —0.5. Perhaps the most impor- 
tant point is that the mean [C/O] ratio in the quasar 
absorption line systems is significantly higher than the 
minimum value seen in metal-poor stars. Unless system- 
atic errors remain in some or all of the measurements, 
this difference may empirically suggest variations in the 
IMF, regardless of theoretical yields. 

Finally, we emphasize that we have compared our re- 
sults mainly to metal-poor stars that are not carbon 
enhanced. These are undcrrepresented in the Cayrel 
et al. (2004) sample, yet their generally poor agreement 
with quasar absorption line abundances is already ap- 
parent from the two such stars included in Figure 12. 
If carbon-enhanced stars fairly reflect their native ISM 



abundances, then these abundances are no longer com- 
mon by z ~ 6. This raises the intriguing possibility that 
most carbon-enhanced stars were formed at even earlier 
times. On the other hand, our high-redshift sample is rel- 
atively small, and we could have easily missed a carbon- 
enhanced population even if it made up ^20% of the to- 
tal number density, similar to the frequency among very 
metal-poor stars (Lucatello et al. 2006). Further surveys 
for high-redshift and metal-poor absorption line systems 
should better establish what fraction of absorbers arc 
carbon-enhanced (e.g., Cooke et al. 2011a), and clarify 
their connection to metal-poor stars. 

6. CONCLUSIONS 

We have measured the relative abundances of C, O, 
Si, and Fe in nine low-ionization quasar absorption line 
systems at 4.7 < z < 6.3. The new data presented here 
are derived from optical and infrared spectra taken with 
Keck, Magellan, and the VLT, and supplement previous 
measurements made with high-resolution optical spectra. 
The absorption lines are generally weak, which allows us 
to measure column densities even for carbon and oxygen, 
which are generally saturated in typical lower-redshift 
systems. 

The relative abundances in our z > 4.7 systems agree 
well with those measured from sub-DLAs and metal-poor 
DLAs at lower redshifts (2 < z < 4.5). The lower- 
redshift sample includes very metal poor systems, where 
carbon and oxygen abundances can be measured (Cooke 
et al. 2011b, and references therein). Among our systems 
and the literature results there is remarkably little appar- 
ent intrinsic scatter (<7q < 0.1 dex) in the column density 
ratio of any two of the dominant ions (C II, O I, Si II, and 
Fe II), and no apparent evolution with redshift. The lack 
of scatter suggests that the relative element abundances 
are correctly estimated from the ratios of the correspond- 
ing ions, with minimal ionization effects or dust deple- 
tion. It also suggests that the metal inventories of most 
metal poor absorption systems over 2 < z < 6 are dom- 
inated by a similar stellar population. At z ~ 6, the fi- 
nite age of the universe means the metals must have been 
produced by short-lived, massive stars. A reasonable in- 
terpretation, therefore, is that the relative abundances 
in metal-poor quasar absorption-line systems reflect the 
integrated yields from prompt supernovae (Type II, Ib/c, 
and prompt la), with minimal contribution from delayed 
Type la supernovae or AGB winds. From these data we 
can directly infer, for example, that massive stars yield 
an iron abundance with respect to oxygen that is roughly 
one third of the solar value. 

The measurements presented here should help in un- 
derstanding the nature of stellar populations at both high 
and low redshifts. The lack of exotic abundance pat- 
terns at z ~ 5-6 suggests that ordinary Pop II (rather 
than massive Pop III) stars are likely to have produced 
most of the ionizing photons during hydrogen rcioniza- 
tion. The broad agreement between the abundances in 
our high-redshift absorption systems and those in metal- 
poor halo stars (e.g., Cayrel et al. 2004) further provides 
a direct link between these stars and the high-redshift 
material out of which they are likely to have formed. 
Future measurements at high redshifts should help to 
strengthen this connection between the oldest stars in 
our galaxy and star formation in the early universe, as 
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well as reveal the metal abundances present at the earli- 
est stages of galactic chemical evolution. 
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TABLE 4 

Column Densities 





QSO 


^abs 


log iV H i 


log N c u 


log N Q 1 


log ^Mg II 


log N Si u 


log N Fe ii 




b h 


Rcf. c 


















(kms" 1 ) 


(km s ) 




SDSS 


J0040-0915 


4.7393 


< 20.4 


> 14.6 


> 15.0 


> 13.4 


14.13 ±0.02 


13.77 ±0.06 


-34,80 




1 


cncc 

OUOO 






< zu.o 




^> I A 7 


i q «c -|_n 11 








Q k _(_ n o 




SDSS 


J0231-0728 


5.3380 d 


< 20.1 


13.79 ±0.05 
14.13 ±0.04 


14.47 ±0.05 


13.02 ±0.01 
13.38 ±0.01 


13.15 ±0.04 
13.42 ± 0.02 


12.73 ±0.04 
12.94 ±0.05 


-25,36 (-38,49) 
-119,160 (-147,184) 




1 


SDSS 


J0818+1722 


5.7911 d 




14.13 ±0.03 
14.19 ±0.03 


14.54 ±0.03 




13.36 ±0.04 
13.41 ±0.03 


12.89 ±0.07 
12.95 ±0.07 


-86,45 (-100,60) 
-124,45 (-140,60) 




1,2 


SDSS 


J0818+1722 


5.8765 




13.62 ± 0.07 


14.04 ±0.05 




12.78 ± 0.05 


< 12.7 


-16,12 (-30,25) 




1.2 


SDSS 


J1148+5251 


6.0115 




14.14 ±0.06 


14.65 ± 0.02 




13.51 ±0.03 


13.21 ±0.25 


-36,32 (-50,45) 




1.2 


SDSS 


J1148+5251 


6.1312 




13.88 ±0.15 


14.79 ± 0.23 




13.29 ±0.06 


13.03 ±0.18 


5.7 ±0.7 


1,2 


SDSS 


J1148+5251 


6.1988 




12.90 ±0.11 


13.49 ±0.13 




12.12 ±0.05 




-12,10 




1.2 


SDSS 


J1148+5251 


6.2575 




13.78 ±0.21 


14.12 ±0.15 




12.92 ± 0.08 e 


< 12.7 




3.8 ±0.3 


1,2 



a Minimum and maximum velocities, with respect to the nominal redshift, over which the optical depths were integrated for systems 
measured using the Apparent Optical Depth method. Where two sets of numbers are given, the first set refers to the HIRES data, and the 
set in parentheses refers to the lower-resolution (X-Shooter or NIRSPEC) data. 

b Dopplcr 6-parameter for systems fit using a Voigt profile 

c References for column density measurements: 1-This work, 2-Becker et al. (2011b) 

d The first line gives the column densities integrated over the velocity interval spanned by O I, while the second line gives the total column 
densities integrated over the entire interval where Si II and C II are detected. 
c Sce notes in text. 
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